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ABSTRACT: In this article, we consider an approach for increasing the density of field-effect transistors, which 
comprised of a digitally controlled oscillator. Based on this approach, we describe manufacture of this oscillator in 
a heterostructure with specialized structure. Doping of some specific regions of the considered heterostructure 
should be done by using diffusion type of doping or by ion type of doping. After doping dopant and radiation 
defects will be annealed by used optimized procedure. We also will consider an approach of decreasing of 
mismatch-induced stress value in this heterostructure. We introduce an analytical approach for modeling of mass 
and heat transport in multilayer structures during manufacturing of integrated circuits taking into account 


mismatch-induced stress. 
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1. INTRODUCTION 

Currently some actual problems of electronics of the 
solid state (one of them is increasing of performance of 
diodes, field-effect and bipolar transistors; another one 
is increasing of their reliability; the third problem is 
increasing of density of elements of integrated circuits) 
are intensively solving [1-7]. One of them is increasing 
of integration rate of elements of integrated circuits, 
because the increasing is one way to decrease their 
dimensions of the considered devices. They are several 
approaches could be used for the increasing. Laser and 
microwave types of annealing of dopants could be used 
to increase the integration rate due to generation of 
inhomogenous distribution of temperature (in this 
situation one can obtain inhomogenous distribution of 
several parameters: diffusion coefficient and other). 
Both types of annealing could be also used to annealed 
radiation defects, generated during ion implantation 
[9-18]. An alternative approach to manufacture of 
integrated circuits with decreased dimensions of 
elements is manufacturing of these elements in this 
film heterostructures by ion implantation or by 
diffusion [1-3]. At the same time it should be noted, 
that the alternative approach requires optimization of 
dopant and/or radiation defects annealing [18]. 


In this paper, we consider an approach of 
manufacturing of field-effect transistors in thin film 
heterostructures, which gives a possibility to decrease 
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dimensions of the above transistors and at the same 
time to increase their density comprised of the digitally 
controlled oscillator [4]. We also described decreasing 
of mismatch-induced stress in the considered 
heterostructure. In this paper framework we will 
consider a heterostructure. The heterostructure consist 
of a substrate and an epitaxial layer (see the 
heterostructure in Fig. 1). We also will consider using 
of a buffer layer between the considered substrate and 
epitaxial layer. By consider the epitaxial layer, which 
includes into itself several sections. These sections 
were manufactured by using other materials. We also 
consider two types of doping of the above sections: by 
dopant diffusion and by ion implantation. One can 
obtain required types of conductivity (p or n) in these 
sections. After that these sections could be used as 
sources, drains and gates (see Fig. 1). Using these ways 
of manufacture of field-effect transistors requires 
optimization of annealing of dopant and/or radiation 
defects. The main aim of this paper is prognosis of 
dopant and radiation defects redistribution to 
determine conditions of decreasing of dimensions of 
considered transistors framework of the oscillator with 
simultaneous increasing of their density. We also 
obtained condition for decreasing of mismatch-induced 
stress in heterostructures. Simultaneous aim of the 
present paper is introduction of an analytical approach 
for prognosis of mass and heat transport in 
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heterostructures during manufacturing of integrated 
circuits with account nonlinearity of considered 
processes and mismatch-induced stress. The analytical 
approach gives us a possibility to analyze mass and 
heat transport in multilayer structures without 
crosslinking of solutions of considered equations on 
interfaces between layers of these structures. The 
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analytical approach also leads to possibility for taking 
into account variation of parameters of considered 
processes in space and time. Using of the approach 
gives us a possibility to analyze manufacturing of the 
considered integrated circuit in more details in 
comparison with recently published works. 


D. O C(x, y,Z,t) is 


Oy Oz OZ 
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Fig. 1. The considered oscillator [4] 
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Function C(x,,z,t) describes the distribution of 
concentration of dopant in space and time; parameter 
Q describes the atomic volume of the considered 
dopant; operator Vs describes symbol of the surficial 


L, 
gradient; integral (C(x, y,Z;, t)d Z describes the 
0 


surfice dopant concentration near interface between 
layers of the considered heterostructure (in this 
situation we assume, that Z-axis is perpendicular to 
interface between layers of heterostructure); function 
ut (xy,z,t) describes the chemical potential, which we 


D, = D,(%,y,z,T) 


Functions Di(x,y,z,T) describes dependence of dopant 
diffusion coefficient on the spatial coordinates (due to 
accounting all layers of heterostruicture) and 
temperature of annealing (due to Arrhenius law); 
parameter T describes annealing temperature; function 
P (x,y,z,T) describes dependence of the limit of 
solubility of dopant the spatial coordinates and 
temperature of annealing; parameter vy depends on 
properties of materials framework considered 
heterostructure and will be integer framework the 
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considering due to the presence of mismatch-induced 
stress in the considered heterostructure; parameters D 
and Ds describe coefficients of volumetric and surface 
diffusions, respectively. Values of the above coefficients 
of diffusions depend on properties of the considered 
heterostructure's materials, materials heating and 
cooling speeds during annealing and distribution of 
dopant concentration in space and time. We 
approximate dependences of considered coefficients of 
diffusions of dopant on considered parameters by the 
following relations (see Refs. [20-22]) 


V"(x,y,z,t) 
(v'y 


following interval y €[1,3] [21]; function V (x,,z,t) 
describes distribution of concentration of radiation 
vacancies in space and time with the equilibrium 
distribution V*. Dependence of diffusion coefficient of 
dopant on dopant's concentration was described in 
details in Ref. [21]. Distributions of concentration of 
point radiation defects in space and time were 
calculated as solution of the following system of 
equations [21,22] 
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Function J (x,y,z,t) describes distribution of 
concentration of radiation interstitials in space and 
time with equilibrium distribution J/*; functions 
Di(x,y,z,T), Dv(xy,z,T) describe dependences of diffusion 
coefficients of interstitials and vacancies on spatial 
coordinates and temperature; quadric terms V2(x,y,z,t) 
and /2(x,y,z,t) of Eqs. (3) gives a possibility to take into 
account generation of divacancies and diinterstitials, 
respectively (see, for example, Ref. [22] and 
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kiv(xy,z,T), ki(xy,z,T) and kvv(xy,z,T) describe 


dependences of parameters of recombination of point 
radiation defects and generation of their complexes on 
spatial coordinates and temperature. 


Distributions of concentrations of divacancies @y 
(xy,z,t) and diinterstitials @® (x,y,z, t) could be 
determined as solutions of the following system of 
equations [21,22] 


appropriate references in this book); functions 
O®,(x, y,z,t) O OD (x, y,z,t) OD (x, yyz,2) 
—Sa ee ee D eae) i rs ==. D 99 a rs 
os Fy Poul 92T J — PS Ay Le) 3 
: Da lo, ,7 Pee) +k, (x, y.2T P(x, y.2.t)—k, (x, ¥,27 Ue, y,25t) 
& z 
(6) 
OD (x, y, zt) O O®D (x, y, zt) O O®D (x, y,z,t) 
. = D an a - DeAXVi2,4 - 
oF os pve) = ee Lene ce 


O OD (x, y,z,t) 
ae D 9N9 wi ¥ 


Boundary and initial conditions for these equations are 


Oo, (x, y,z,t) 


Kgs (x, y, z,T V"(x, y, z,t)— k, (x, y, z,T)V (x, y, z,t) 


Ox Ox 


x=0 


0® | (x, y,z,t) 


dy 


O® , (x, y,z,t) 
OZ 


, 


C2 


0. O® (x, y,z,t) 


_0 OO, (x, y,z,t) 


O® (x, y,Z,t) 
: By 


=0 


x=Ly 


=(0, 


y=0 


=) 


= (0), D(x,y,z,0)=for (x,y,z), Dr (x,y,z,0)=fov(xy,z). (7) 


Pankratov 


International Journal of Advanced Science and Engineering 


www.moahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.7 No.2 1652-1674 (2020) 1656 


Here Dai(x,y,z,T) and Dov(x,y,z,T) are the coefficients of 
diffusions of complexes of radiation defects; ki(x,y,z,T) 
and ky(x,y,z,T) are the parameters of decay of 
complexes of radiation defects. 
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point radiation defects, 0<¢<1, |g,(xy,z,T)|<1, e =LV. 
One can also use analogous transformations of 
approximations of recombination parameters of point 
radiation defects and parameters of their complexes 


generation: kiv(X,y,Z,T)=Korv[1+é,v givlXy,z,T)], 
ki(xy,Z,T)=ko1 [1+ 61 gi(Xy,z,T)| and kv v(x,y,Zz,T)=kovv 
[1+evv gvv(xy,z,T)]/ Here kopi,2 describe average 
values of the considered, 0<é,v<1, O0<é <1, O<evv<1, | 


giWX&y,Z,T)|<1, | gulxyz,T)|<1, |gvv(xy,z,T)|<1. Let us 
introduce the following dimensionless variables: 
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Now we will calculate solutions of the obtained Eqs.(8) with conditions (9) by using recently introduced approach 
[18], i.e. we consider the required solution as the following power series 
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P(yn.g,9)= Ye, LO! LOD, (7.7.6.9), (10) 


Now let us substitute the above series in the Eqs. (8) 
and (9). The substitution leads to obtain equations for 


initial-order approximations of concentration of 
considered defects I ( 2: 1, ?, 9) and 
Vino XM ?, 9). At the same time we obtain equations 
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i 21, j 21, k2=1. We present these equations in Appendix 
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standard approach of Fourier [24,25]. We presented all 
obtained solutions in Appendix for the paper. 


Next we will calculate spatio-temporal distributions of 
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concentrations. First of all let us transform 
approximations of diffusion coefficients in the 
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Here Dowp are the average values of the considered 
diffusion coefficients. By substitution of transformation 
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Let us to substitute series (11) into Eqs.(6) and 
appropriate conditions: both boundary and initial. The 
substitution leads to obtaining equations for 
calculation the initial-order approximations of simplest 
complexes of radiation defects concentrations 
@p»(x,y,z,t), all corrections for the _ initial-order 
approximations @,(x,y,z,t) (for them i 21) and 
conditions for them (both: boundary and initial). We 
remove the obtained equations and conditions to 
Appendix for the paper. Solutions of the considered 
equations were calculated by using standard Fourier 


approach [24,25] and presented in the considered 
Appendix in end of the paper. 


Next let us calculate dopant concentration as a 
distribution in space and time. We used again the same 
approach, which we used for analysis of distributions 
of concentrations of radiation defects. Based on the 
approach we use the following transformation of 
approximation of diffusion coefficient of dopant: 
Di(X,y,Z,T)=Do1[1+ 21g1(Xy,Z,T)]. Parameter Do. 
describes the average value of diffusion coefficient of 
dopant, O<a< 1, |gi(xy,z,T)|<1. After using the 
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form 
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Now let us substitute the considered series in Eq.(1) 
and conditions (2). The substitution gives a possibility 
to obtain equations for functions Cj(x,,z,t) at any 
values of indexes i and j, initial and boundary 
conditions for these functions. The obtained equations 
were presented in the Appendix in the end of this 
paper. We solved the obtained equations by standard 
Fourier approach (see Ref. [24,25]). The obtained 
solutions were presented in the same Appendix in the 
end of this paper. 


Framework this paper we make analysis of 
distributions of dopant and_ radiation defects 
concentrations. These distributions were considered as 
a function of spatial coordinate, current time and other 
parameters. The analysis has been done analytically: 
we used the second-order approximations on all 
parameters, which have been used in appropriate 
series. Usually the second-order approximations are 
enough good approximations to make qualitative 
analysis and to obtain quantitative results. All 
analytical results have been checked by numerical 
simulation. 


DISCUSSION 


Framework the section we make analysis of the spatio- 
temporal distributions of dopant concentration with 
account redistribution of radiation defects. We present 
on Figs. 2 some typical distributions in space of dopant 
concentrations near interfaces of the considered 
heterostructures. We _ obtain’ the considered 
distributions at the condition: diffusion coefficient of 
dopant in doped area has larger value in comparison 
with analogous values in the nearest areas. These 
figures show that in this case one can obtain decreasing 
of dimensions of the considered transistors. At the 
same time one can find increasing of homogeneity of 
distribution of dopant concentration. One can change 
relation between values of dopant diffusion 
coefficients. But this changing leads to opposite result. 
Results of this changing were presented on Figs. 3. 
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perpendicular to interface between materials of 
heterostructure. The considered distributions were 
calculated at the condition: diffusion coefficient of 
dopant in doped area has larger value in comparison 
with analogous values in the nearest areas. Increasing 
of number of curves corresponds to increasing of 
difference between values of diffusion coefficient of 
dopant in layers of heterostructure 
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Fig 2b. Spatial distributions of implanted dopant 


concentration. Direction of the distribution is 
perpendicular to interface between materials of 
heterostructure. Increasing of number of curves 


corresponds to increasing of difference between values 
of diffusion coefficient of dopant in layers of 
heterostructure. Curve 1 corresponds to homogenous 
sample and annealing time © = 0.0048 (Ly?+L,?+ 
L*)/Do. Curve 2 corresponds to homogenous sample 
and annealing time © = 0.0057 (L,?+L,?+ Lz*)/Do. 
Curves 3 and 4 correspond to heterostructure from Fig. 
1; annealing times © = 0.0048 (L,?+L,*+L2*)/Do and © = 
0.0057 (L;?+Ly?+Lz*)/Do, respectively 
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Fig. 3a. Spatial distributions of infused dopant 
concentration in average section of epitaxial layer of 
the considered heterostructure. The dopant was 


0 L/4 L/2 3L/4 I infused in direction, which is parallel to interface 
’ between epitaxial layer and _ substrate’ of 
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heterostructure. Increasing of number of curves 
corresponds to increasing of difference between values 
of diffusion coefficient of dopant in layers of 
heterostructure. Value of dopant diffusion coefficient in 
this section is smaller, than value of dopant diffusion 
coefficient in nearest sections 
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Fig.3b. Spatial distributions of implanted dopant 
concentration in average section of epitaxial layer of 
the considered heterostructure. The dopant was 
infused in direction, which is parallel to interface 
between’ epitaxial layer and_ substrate’ of 
heterostructure. Increasing of number of curves 
corresponds to increasing of difference between values 
of diffusion coefficient of dopant in layers of 
heterostructure. Value of dopant diffusion coefficient in 
this section is smaller, than value of dopant diffusion 
coefficient in nearest sections. Dushed lines are spatial 
distributions of implanted dopant in one epitaxial layer 
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Fig.4a. Distributions of concentration of infused dopant 
in depth of heterostructure from Fig. 1 for different 
values of annealing time (curves 2-4) and idealized 
step-wise approximation (curve 1). Increasing of 
number of curve corresponds to increasing of 
annealing time. 
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Fig. 4b. Distributions of concentration of implanted 
dopant in depth of heterostructure from Fig. 1 for 
different values of annealing time (curves 2-4) and 
idealized step-wise approximation (curve 1). 
Increasing of number of curve corresponds to 
increasing of annealing time 
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Fig. 5a. Dimensionless optimal annealing time of 
infused dopant as a function of several parameters. 
Curve 1 describes the dependence of the annealing 
time on the relation a/L and €= y = 0 for equal to each 
other values of dopant diffusion coefficient in all parts 
of heterostructure. Curve 2 describes the dependence 
of the annealing time on value of parameter «é for 
a/L=1/2 and € = y = 0. Curve 3 describes the 
dependence of the annealing time on value of 
parameter € for a/L=1/2 and ¢ = y = 0. Curve 4 
describes the dependence of the annealing time on 
value of parameter yfor a/L=1/2 and ¢ = €=0 
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case (after finishing of the annealing) dopant should 
achieve interface between layers of the considered 
heterostructure. If the dopant cannot achieve the 
interface during the annealing, then it is attracted an 
interest to make additional annealing. The Fig. 5b 
shows the same dependences of the considered 
optimal values of additional annealing time as on Fig. 
5a. Necessity to anneal radiation defects leads to use 
shorter optimal annealing of implanted dopant in 
comparison with same annealing of the infused dopant. 
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Fig. 5a. Dimensionless optimal annealing time of 
implanted dopant as a function of several parameters. 
Curve 1 describes the dependence of the annealing 
time on the relation a/L and = y = 0 for equal to each 
other values of dopant diffusion coefficient in all parts 
of heterostructure. Curve 2 describes the dependence 
of the annealing time on value of parameter ¢ for 
a/L=1/2 and € = y = 0. Curve 3 describes the 
dependence of the annealing time on value of 
parameter € for a/L=1/2 and ¢ = y = 0. Curve 4 
describes the dependence of the annealing time on 
value of parameter yfor a/L=1/2 and ¢ = €=0 
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Fig. 5a. Dimensionless optimal annealing time of 
implanted dopant as a function of several parameters. 
Curve 1 describes the dependence of the annealing 
time on the relation a/L and €= y = 0 for equal to each 
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other values of dopant diffusion coefficient in all parts 
of heterostructure. Curve 2 describes the dependence 
of the annealing time on value of parameter ¢ for 
a/L=1/2 and € = y = 0. Curve 3 describes the 
dependence of the annealing time on value of 
parameter € for a/L=1/2 and ¢ = y = 0. Curve 4 
describes the dependence of the annealing time on 
value of parameter yfor a/L=1/2 and ¢ = €=0 


CONCLUSIONS 


In this paper we introduce an approach to increase 
integration rate of element of a digitally controlled 
oscillator. The approach gives us possibility to decrease 
area of the elements with smaller increasing of the 
element’s thickness. 
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Solutions of the above equations, obtained by standard Fourier, could be written as 
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